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Abstract 

The  interaction  realized  by  fuel  cell — microturbine  hybrids  derive  primarily  from  using  the  rejected  thermal  energy  and  combustion  of 
residual  fuel  from  a  fuel  cell  in  driving  the  gas  turbine.  This  leveraging  of  thermal  energy  makes  the  high  temperature  molten  carbonate  fuel 
cells  (MCFCs)  ideal  candidates  for  hybrid  systems.  Use  of  a  recuperator  contributes  to  thermal  efficiency  by  transferring  heat  from  the  gas 
turbine  exhaust  to  the  fuel  and  air  used  in  the  system. 

Traditional  control  design  approaches,  consider  a  fixed  operating  point  in  the  hope  that  the  resulting  controller  is  robust  enough  to  stabilize 
the  system  for  different  operating  conditions.  On  the  other  hand,  adaptive  control  incorporates  the  time-varying  dynamical  properties  of  the 
model  (a  new  value  of  gas  composition)  and  considers  the  disturbances  acting  at  the  plant  (load  power  variation). 

©  2002  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

Fuel  cells  are  particularly  well  suited  for  power  generation 
in  co-generation  plants  because  they  convert  energy  directly 
into  electricity  in  an  electrochemical  process  while  simul¬ 
taneously  producing  heat  [1-3].  However,  it  is  first  neces¬ 
sary  to  reform  fuels  such  as  natural  gas,  converting  it  into  a 
gas  with  a  high  hydrogen  content,  so  that  it  can  be  electro- 
catalytic  ally  oxidized  with  air. 

The  direct  carbonate  fuel  cell  is  a  variant  of  molten 
carbonate  fuel  cells  (MCFCs)  in  that  it  internally  reforms 
methane-containing  fuels  within  the  anode  compartment  of 
the  fuel  cell  [4-7],  The  largest  demonstration  of  MCFC 
technology  has  been  California’s  2-MW  Santa  Clara 
Demonstration  Project  [8]. 

Microturbines,  which  are  typically  fueled  with  natural 
gas,  generate  between  25  and  200  kW  of  electricity.  Their 
small  size  and  relatively  low  cost  allow  them  to  be  located 
near  where  they  are  needed.  They  can  operate  at  very  low 
emission  levels  and  reduce  the  efficiency  losses  and  envir¬ 
onmental  impact  of  large  transmission  and  distribution 
systems.  In  this  paper,  a  MCFC  is  associated  with  a  gas 
microturbine  to  produce  electric  power. 
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Hybrid  systems  offer  a  solution  to  two  important  pro¬ 
blems,  the  low  efficiency  and  relatively  high  emissions  of 
small  gas  turbines,  and  the  high  cost  of  small  fuel  cell  power 
plants  [9-11]. 

The  paper  is  structured  as  follows.  Section  2  presents  a 
review  of  the  MCFC.  Some  basic  concepts  of  the  gas  turbine 
theory  are  presented  in  Section  3.  Section  4  describes  the  gas 
turbine  control  configuration.  Section  5  briefly  discusses  the 
fuel  cell — microturbine  hybrid  power  cycles.  Section  6 
outlines  the  adaptive  control.  Section  7  depicts  some  simu¬ 
lation  results  and  discussion.  Finally,  conclusions  are  pre¬ 
sented  in  Section  8. 


2.  Molten  carbonate  fuel  cells 

MCFC  can  reach  fuel-to-electricity  efficiencies  approach¬ 
ing  60%,  considerably  higher  than  the  37-42%  efficiencies 
of  a  phosphoric  acid  fuel  cell  (PAFC)  plant.  When  the  waste 
heat  is  captured  and  used,  overall  fuel  efficiencies  can  be  as 
high  as  85%. 

Improved  efficiencies  are  one  reason  why  MCFC  offers 
significant  cost  reductions  over  PAFC  technology.  Another 
is  that  the  electrodes  of  a  MCFC  be  made  of  nickel  catalysts 
rather  than  the  more  costly  platinum  of  PAFC  systems. 

Natural  gas  is  internally  reformed,  partially  in  an  internal 
reforming  unit  and  partially  at  the  cells,  eliminating  the  need 
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Nomenclature 

Fuel  cell 

E0 

standard  potential 

F 

Faraday’s  constant 

i 

cell  load  current 

Pi 

partial  pressure,  Pia  =  xaiPa/Pa tm, 

P  i,  c  =  XciPc/ P  atm 

Vc 

cell  voltage  under  load  current 

Vo 

equilibrium  potential 

z 

cell  ohmic  impedance 

fad 

activation  polarization 

fconc 

concentration  polarization 

Gas  turbine 

a ,  b,  c 

valve  parameters 

Fpa 

specific  heat  of  air  at  constant  pressure 

CPg 

specific  heat  of  combustion  gases 

Cps 

specific  heat  of  steam 

«1 

valve  position 

Fd 

fuel  demand  signal 

Ahic 

isentropic  enthalpy  change  for  a  compres¬ 
sion  from  pc.m  to  pCout 

A  hyY 

isentropic  enthalpy  change  for  a  gas 
expansion  from  pT|n  to  pT()ut 

Ah25 

specific  enthalpy  of  reaction  at  reference 
temperature  of  25  °C 

HHV 

higher  heating  value 

kf 

fuel  system  gain  constant 

^HHV 

factor  which  depends  on  HHV 

Pcin 

air  pressure  at  compressor  inlet 

Pc  out 

air  pressure  at  compressor  outlet 

PTm 

pressure  of  combustion  gases  at  turbine  inlet 

PTout 

pressure  of  combustion  gases  at  turbine 
outlet 

Pc 

compressor  power  consumption 

p 

1  mec 

mechanical  power  delivered  by  turbine 

Pr 

total  mechanical  power  delivered  by 
turbine 

t 

time 

T 

mechanical  torque  delivered  by  turbine 

T’cou, 

outlet  air  temperature 

T1S 

temperature  of  injected  steam 

Ft 

1  in 

turbine  inlet  gas  temperature 

wa 

air  mass  flow  into  the  compressor 

Wf 

fuel  mass  flow 

Wg 

turbine  gas  mass  flow 

Wis 

injection  steam  mass  flow 

fc 

overall  compressor  efficiency 

f  trans 

transmission  efficiency  from  turbine  to 
compressor 

fT 

overall  turbine  efficiency 

Tf 

fuel  system  time  constant 

CO 

rotation  speed  of  the  turbine 

Aco 

rotation  speed  deviation  of  the  turbine 

Adaptive  control 

A(z~),  B(z~) 

polynomials  of  the  plant 

polynomials  of  the  performance  index 

B~(z~) 

e(i) 

disturbance  of  the  plant 

h 

constant  term 

Aide 

stack  current  deviation 

J 

performance  index 

K 

prediction  horizon  in  the  performance 
index 

filter  polynomials 

Viz-1) 

P  init 

covariance  matrix 

apl 

load  power  variation 

q 

weighting  coefficient 

r\,  r2 

varying  parameters  of  the  polynomial 
Pcizr') 

Riz~\  S{zT\ 

polynomials  of  the  difference  equation 

Tiz-1) 

Vi 

standard  deviation 

51>  s2 

varying  parameters  of  the  polynomial 
Sc(z-‘) 

Sciz~\  Rc(z~l) 

polynomials  of  the  controller 

u(i) 

plant  control  signal 

yd ) 

plant  output  signal 

2 

forgetting  factor 

(7~ 

constant  variance 

a)(i) 

setpoint  of  the  plant  output  signal 

for  a  large  external  reforming  unit  to  produce  hydrogen  fuel. 
The  approach  (Fig.  1)  is  a  combination  of  indirect  internal 
reforming  (HR)  and  direct  internal  reforming  (DIR)  which 
provides  for  better  thermal  management.  At  the  cathode, 
oxygen  reacts  with  carbon  dioxide  and  electrons  to  form 
carbonate  ions: 

^02  +  C02  +  2e“  — >  C032-  (1) 

The  carbonate  ions  flow  through  the  electrolyte  matrix 
from  cathode  to  anode.  At  the  anode,  the  carbonate  ions  are 
consumed  by  the  oxidation  of  hydrogen  to  form  steam  and 
carbon  dioxide,  releasing  electrons  to  the  external  circuit: 

H2  +  C032-  -»•  H20  +  C02  +  2e"  (2) 

Cell  voltage  under  load  current  (A/cm)  is  expressed  as 

[12]: 

Vc  V)  17act  f  cone  E-  (^) 

Activation  polarization  is  caused  by  electrode  kinetics 
while  concentration  polarization  is  caused  by  concentration 
gradients  in  the  electrode.  Equilibrium  potential  is  described 
by  the  Nernst  equation  [12]: 


V0  =  E0  +  —  in 
2  F 


Pn20,aPc02,a 


(4) 
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Fig.  1.  IIR/DIR  structure  of  MCFC  stack. 


Partial  pressures  (normalized  to  atmospheric  pressure) 
depend  on  anode/cathode  gas  pressure  and  composition  while 
standard  potential  and  ohmic  impedance  are  both  tempera¬ 
ture  dependent.  Fuel  cell  polarization  losses  are  generally 
dependent  on  partial  pressures,  temperature  and  current 
density,  and  are  spatially  distributed  in  an  actual  cell.  In  this 
paper,  the  dynamic  model  is  a  lumped  parameter,  where  outlet 
properties  are  equal  to  average  properties  [13-15]. 


3.  Gas  turbine 

A  thorough  introduction  to  gas  turbine  theory  is  provided 
in  the  book  of  Cohen  et  al.  [16].  There  also  exist  a  large 
literature  on  the  modeling  of  gas  turbines.  Model  complexity 
varies  according  to  the  intended  application.  Detailed  first 
principle  modeling  based  upon  fundamental  mass,  momen¬ 
tum  and  energy  balances  is  reported  by  Fawke  et  al.  [17]  and 
Shobeiri  [18].  These  models  describe  the  spatially  distrib¬ 
uted  nature  of  the  gas  flow  dynamics  by  dividing  the  gas 
turbine  into  a  number  of  sections.  Throughout  each  section, 
the  thermodynamic  state  is  assumed  to  be  constant  with 
respect  to  the  location,  but  varying  with  respect  to  time. 
Mathematically,  the  full  partial  differential  equation  model 
description  is  reduced  to  a  set  of  ordinary  differential 
equations,  which  facilitate  easier  application  within  a  com¬ 
puter  simulation  program.  For  a  detailed  model,  a  section 
might  consist  of  a  single  compressor  or  turbine  stage.  Much 
simpler  models  result  if  the  gas  turbine  is  decomposed  into 
just  three  sections  corresponding  to  the  main  turbine  com¬ 
ponents,  i.e.  compressor,  combustor  and  turbine  [19]. 

Another  modeling  approach  is  by  utilizing  real  steady- 
state  engine  performance  data  [20].  It  is  assumed  that 


transient  thermodynamic  and  flow  processes  are  character¬ 
ized  by  a  continuous  progression  along  the  steady-state 
performance  curves,  this  is  known  as  the  quasi-static 
assumption.  The  dynamics  of  the  gas  turbine,  e.g.  combus¬ 
tion  delay,  motor  inertia,  fuel  pump  lag,  etc.  are  then 
represented  as  lumped  quantities  separate  from  the 
steady-state  performance  curves.  Very  simple  models  result 
if  it  is  further  assumed  that  the  gas  turbine  is  operated  at  all 
times  close  to  rated  speed  [21], 

In  terms  of  energy  conversion,  chemical  energy  present  in 
the  combustion  reactants  is  transferred  to  the  gas  stream 
during  combustion.  This  energy — measured  in  terms  of  gas 
enthalpy  or  higher  heating  value  (HHV) — is  then  converted 
into  mechanical  work,  by  expanding  the  gas  through  the 
turbine.  Thus,  the  excess  mechanical  power  available  for 
application  elsewhere,  after  accounting  for  the  power 
required  to  drive  the  compressor  is  derived  ultimately  from 
the  combustion  process.  Without  combustion,  assuming 
100%  efficient  compressor  and  turbine  operation,  the  power 
developed  by  the  turbine  would  be  exactly  matched  by  the 
power  required  to  drive  the  compressor. 

Compressor  power  consumption  equation: 


Pc 


WaA/ric 

^/c^/trans 


(5) 


Combustion  energy  equation: 


wgCpg(rTin  -  298)  +  WfA/225  +  wacPa(298  -  TCout) 

+  wisCpS(298  -  7*)  =  0  (6) 

Power  delivery  equation: 

pT  =  iiTWgAhn  (7) 

Pr nec  =  Pj  Pc  (8) 


Fig.  2  shows  the  block  diagram  of  the  gas  turbine. 


4.  Gas  turbine  control  configuration 

The  concept  of  the  gas  turbine  control  system,  which  is 
applied  in  this  paper,  is  based  on  the  Speedtronic  Mark  4 
description  as  presented  in  the  paper  by  Rowen  [22].  Some 
considerations  concerning  the  subject  may  be  also  found  in 
[23,24], 

The  simplified  gas  turbine  model  is  divided  into  two 
interconnected  subsystems  in  this  paper.  The  subsystems 
are:  the  fuel  system  (fuel  valve  with  actuator),  and  the 
turbine. 

The  fuel  system  consists  of  the  fuel  valve  and  the 
actuator.  The  fuel  flow  out  from  the  fuel  systems  result 
from  the  inertia  of  the  fuel  system  actuator  and  of  the  valve 
positioner. 

Fuel  system  actuator  equation: 


W  f 


k{ 

TfS  +  1 


e\ 


(9) 
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Fig.  2.  Gas  turbine. 


Valve  positioner  equation: 

ci  =t^~Fa  (10) 

Input  variable  to  the  fuel  system  is  Fd.  Output  variable  from 
the  fuel  system  model  is  Wf. 

A  single  gas  turbine  does  not  require  the  digital  setpoint 
feature.  The  turbine  torque  function  is  given  by: 

T  =  ^HHv(wf  —  0.23)  +  0.5  (Am)  (11) 

The  A:Hhv  and  0.23  factors  cater  for  the  typical  turbine 
power/fuel  rate  characteristic,  which  rises  linearly  from  zero 
power  at  23%  fuel  rate  to  rated  output  at  100%  fuel  rate. 

The  Eq.  (11)  allows  the  turbine  torque  to  be  calculated 
algebraically.  This  torque  is  used  in  the  equations  which 
model  the  mechanical  system: 

P  mec  =  Tea  (12) 

Input  variables  to  the  turbine  are  Wf,  Aco  and  co.  Output 
variable  from  the  turbine  is  Pmec. 

For  the  purpose  of  this  paper,  only  modulating  control  of 
mechanical  side  of  the  gas  turbine  is  of  interest.  The 


simplified  model  of  the  gas  turbine  controller  in  this  paper 
consists  of  two  inputs  and  one  output.  Inputs  to  the  controller 
are  Pmec  and  co.  Output  from  the  controllers  is  Fd. 

The  block  diagram  of  the  gas  turbine  control  system  is 
presented  in  Fig.  3.  The  diagram  consists  of  two  PID 
controllers.  FVG  stands  for  least  value  gate  that  transmits 
the  minimum  of  two  incoming  signals. 

5.  Fuel  cell/turbine  hybrid 

When  used  in  combination  with  turbines,  fuel  cells  can 
produce  from  55  to  90%  of  the  electricity  of  the  system 
while  turbines  produce  the  remainder  [10].  Several  cycle 
configurations  have  been  proposed  and  the  terminology  is 
still  evolving,  but  one  useful  way  of  looking  at  the  differ¬ 
ences  is  to  divide  the  configurations  into  those  that  include 
“directly  fired  turbines”  and  those  with  “indirectly  fired 
turbines”.  A  hybrid  power  system  with  a  “directly  fired 
turbine”  normally  uses  a  pressurized  fuel  cell  to  provide 
input  to  the  turbine,  thus,  acting  as  a  combustor.  In  the 
“indirectly  fired  turbine”  system  an  atmospheric  fuel  cell  is 
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co„ 


co 
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Fig.  3.  Gas  turbine  control  system. 
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Fig.  4.  Fuel  cell/turbine  hybrid  system. 


used  and  a  pressurized  heat  exchanger  provides  input  to  the 
turbine.  When  utilizing  the  fuel  cell/gas  turbine  combination 
the  combined  efficiency  of  the  system  is  raised  to  >60%  and 
criteria  pollutant  emissions  are  essentially  eliminated. 

An  example  of  a  hybrid  system  with  an  “indirectly  fired 
turbine”  is  shown  in  Fig.  4.  This  type  of  integration  could 
utilize  any  gas  turbine  system.  The  indirectly  fired  approach 
posits  the  fuel  cell  in  the  exhaust  of  an  indirectly  fired  gas 
turbine.  The  air  and  residual  products  of  the  fuel  cell  are  then 
fed  to  an  atmospheric  combustor.  This  combustor  heats  the 
air  leaving  the  compressor  via  a  heat  recovery  unit  (HRU) 
and  delivers  it  to  the  turbine.  The  combination  of  the 
atmospheric  combustor  and  heat  exchanger  replaces  the 
normal  internal  pressurized  combustor. 

The  turbine  exhaust  flows  to  the  fuel  cell  anode  exhaust 
oxidizer.  Exhaust  from  the  anode  exhaust  oxidizer  flows  to 
the  heat  exchanger,  which  provides  the  heat  for  the  com¬ 
pressor  air.  The  exit  from  the  heat  exchanger  flows  through 
the  fuel  cell  cathode  providing  the  oxygen  and  carbon 
dioxide  needed  in  the  carbonate  fuel  cell  process. 

In  the  “indirectly  fired  turbine”  system,  the  fuel  cell  does 
not  need  to  operate  at  the  turbine  pressure,  instead  it  operates 
at  the  preferred  ambient  pressure  and  is  independent  of  gas 
turbine  cycle  pressure  ratio.  The  system  works  efficiently 
with  a  wide  range  of  turbine  compression  ratio.  This  allows 
taking  a  system  developed  for  integration  at  the  multi-MW 
scale  with  industrial  size  turbines,  and  configuring  a  small 
MW  class  system  using  a  microturbine  at  a  lower  pressure 
ratio. 


6.  Adaptive  control 

A  real-world  plant  can  be  usually  characterized  by 
time-varying  dynamical  properties,  which  affect  the  plant 
behavior.  Stochastic  models  are  used  to  represent  the 
disturbances  acting  at  the  plant  output  because  of  the 


large  number  and  different  nature  of  the  factors  disturbing 
the  normal  plant  operation. 

Robustness  properties  can  usually  be  ensured  by  the 
feedback  structure  of  the  control  system.  The  feedback 
compensates  for  the  deviation  of  the  plant  output  signal 
value  from  its  setpoint:  disturbances  affecting  the  plant  (load 
power  variation)  or  change  in  the  plant  model  parameters 
(HHV),  such  a  change,  is  convenient  to  observe  parameter 
tracking. 

It  may  be  possible  to  identify  the  parameters  of  the 
controller  that  we  are  seeking.  This  scheme  is  called  direct 
adaptive  control,  because  we  are  going  to  obtain  directly  the 
required  controller  parameters  through  their  estimation  in  an 
appropriately  redefined  plant  model. 

Adaptive  control  is  usually  used  to  cope  with  an  unknown 
or/and  changing  plant  to  be  controlled  [25].  Analysis  and 
synthesis  of  such  a  control  system  is  possible  only  under 
some  assumptions  concerning  the  nature  of  the  plant  and  its 
dynamics.  In  this  paper,  only  linear,  discrete-time  plants 
disturbed  in  a  stochastic  manner  will  be  considered.  The 
following  plant  model  will  be  used  [26]: 

y(i)  =  +  e(0  (13) 

Eq.  (13)  is  one  of  the  most  typical  in  the  field  of  adaptive 
control  and  non-standard  discrete-time  control  algorithms 
in  general.  If  only  stable  factors  exist  in  the  B  polynomial, 
the  plant  will  be  called  minimum  phase.  The  part  e(i)  is 
the  stochastic  part  of  the  disturbance.  In  this  paper,  this 
disturbance  is  the  load  power  variation  A PL. 

Generally,  most  control  algorithms  can  be  described  by 
the  structure  and  parameters  of  the  difference  equation  [26]: 

R(z~l)u(i)  +  S(z~l)y(i)  -  T(z~l)oo(i)  +  h  =  0  (14) 

The  coefficients  of  the  Rlz^1),  5(z^!)  and  T(z_1)  poly¬ 
nomials  and  the  h  term  are  chosen  before  the  simulation 
experiment  and  stay  constant  during  the  experiment. 
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Fig.  5.  Fuel  cell  control  system. 


The  aim  of  the  minimum  variance  control  algorithm  is  the 
minimization  of  the  following  performance  index  [25,27]: 


J  =  E<  P(z~')  +  -  v(z->( 0 


+  cj[Q(z  l)u(iY] 


(15) 


where  B,  _(z_1)  is  the  reciprocal  polynomial  of  the  B  (z  ') 
polynomial.  Minimization  of  this  performance  index  leads 
to  a  control  algorithm  of  the  same  structure  as  Eq.  (14). 

The  estimation  scheme  used  in  this  paper  is  the  recursive 
least  squares  (RLS).  The  adaptive  control  is  shown  in  Fig.  5. 
Stack  current  Idc  is  measured  and  used  in  calculation  of 
fuel  flow  setpoint.  This  measurement  of  the  stack  current 
determines  the  plant  output  signal  y(i).  Stack  current  is 
continuously  adjusted  by  the  inverter  control  to  maintain 
power.  This  current  is  the  setpoint  of  the  plant  output  signal 
co(i). 


7.  Results  and  discussion 


jected  to  a  sudden  valuation  in  power  demand.  The  HRU, 
power  conditioning  system  (PCS),  and  plant  control  system 
are  included  in  the  simulation  [28-31].  The  inverter  is 
assumed  to  regulate  load  voltage  perfectly,  and  simply  draws 
stack  current  proportional  to  load  current  and  inversely 
proportional  to  stack  voltage  using  a  power  demand  setpoint. 
The  time  constants  for  changes  in  power  output  for  the 
microturbines  and  fuel  cell  range  from  5  ms  to  50  s,  and  the 
PCS  dynamics  are  not  important  [32], 

This  paper  develops  the  control  system  with  an  adaptive 
minimum  variance  controller.  The  plant  and  the  controller 
are  simulated  as  discrete  in  time  [33].  The  plant  parameters 
are:  polynomials  A(z~  '),  B(z  ')  and  standard  deviation  ,vd  of 
white  noise  of  normal  distribution.  Parameters  of  RLS  are: 
initial  values  of  diagonal  of  covariance  matrix  P_in it  and 
forgetting  factor  L 

The  disturbance  considered  is  a  load  power  variation  A PL. 
This  is  zero-mean,  white  noise  and  has  constant  variance  a  . 

The  controller  has  the  form  [27]: 

“(0  =  _IH3’(i)  (16) 


A  plant  consisting  of  a  load  is  fed  from  the  fuel  cell/ 
turbine  hybrid  system.  The  selected  system  comprises  a 
250  kW  fuel  cell  and  a  30  kW  gas  microturbine.  The  plant 
and  the  fuel  cell/microturbine  system  are  modeled  using 
MATLAB™. 

All  parameters  correspond  to  a  two-stack  equivalent.  The 
fuel  cell  stacks  used  in  this  paper  are  rated  at  125  kW.  Stack 
voltage  is  taken  across  a  parallel  connection  of  two  stacks, 
each  stack  consisting  of  258  cells. 

To  investigate  transient  behavior,  the  plant  is  assumed  to 
be  at  steady-state  corresponding  to  rated  power  and  sub¬ 


where  Sc  and  Rc  are  of  order  2.  The  structure  of  the  controller 
is  calculated  according  to  the  plant  structure.  To  show  the 
ability  of  the  controller  to  adapt  to  varying  operating  con¬ 
ditions,  a  new  value  of  gas  composition  is  introduced  at 
100  s  what  provides  time-varying  plant.  Controller  para¬ 
meters  are  identified  on-line  and  the  simulation  time  is  200  s. 

The  minimum  variance  controller  serves  as  an  example  of 
a  self-tuner.  The  optimal  value  of  the  output  signal  standard 
deviation  is  obtained  in  the  steady-state.  The  plant  is 
assumed  to  be  of  Auto  Regressive  with  exogenous  input 
(ARX)  type. 


Current  (A)  S  ln,Put  sequence 
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Fig.  6.  (a)  Control  signal  U\  (b)  load  current  deviation  A/dc. 
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The  results  with  the  derived  model  are  summarized. 
Fig.  6a  presents  the  control  signal  U.  Referring  to  the  results 
shown  in  Fig.  6b,  it  can  be  seen  the  stack  current  deviation 
A/dc.  Stack  voltage  starts  at  its  initial  value,  corresponding  to 
regulated  power,  and  drops  quickly  to  a  sudden  increase  in 
stack  current,  in  turn,  controlled  by  the  inverter  to  maintain 
the  new  power  setpoint.  Stack  voltage  varies  according  to  the 
temperature  dependence,  while  stack  current  is  continuously 
adjusted  by  the  inverter  control  to  maintain  power. 

Fig.  7  shows  the  output  standard  deviation  ,vd  and  the 
optimal  output  standard  deviation.  Fig.  8a  depicts  the  poly¬ 
nomial  R  parameters  where  the  varying  parameters  consid¬ 
ered  are  rx  and  r2.  Fig.  8b  displays  the  polynomial  S 
parameters,  .v,  and  s2  being  the  varying  parameters.  This 
controller  guarantees  some  additional  robustness  margins  in 
the  case  that  the  model  does  not  cover  the  entire  plant 
uncertainty. 

8.  Conclusion 

The  emerging  fuel  cells  produce  very  high  temperature 
exhaust  gases  that  can  be  used  to  drive  a  gas  turbine.  A 
dynamic  model  for  this  MCFC — microturbine  system  has 
been  elaborated. 

A  MCFC  stack  interacts  with  other  system  components 
and  causes  an  interdependency  between  them.  To  investigate 
transient  behavior,  the  plant  is  assumed  to  be  at  steady-state 
corresponding  to  rated  power  and  subjected  to  a  sudden 


variation  in  power  demand.  This  paper  develops  the  control 
system  with  an  adaptive  minimum  variance  controller  and 
based  on  the  simulation  study,  the  resulting  controller  is 
robust  enough  to  stabilize  the  system  for  different  operating 
conditions. 
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